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Abstract 

We propose a model based on quantum molecular dynamics (QMD) incorpo- 



rated with statistical decay model (SDM) to describe various nuclear reactions 



^ ' in an unified way. In this first part of the work, the basic ingredients of the 



model are defined and the model is applied systematically to the nucleon(iV)- 
induced reactions. It has been found that our model can give a remarkable 
agreement in the energy-angle double differential cross sections of (TV, xN 1 ) 
type reactions for incident energies from 100 MeV to 3 GeV with a fixed pa- 
rameter set. An unified description of the major three reaction mechanisms of 
(N,xN') reactions, i.e. compound, pre-equilibrium and spallation processes, 

is given with our model. 
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I. INTRODUCTION 



Nuclear reactions reveal various aspects of hadronic many-body problem as a function 
of the target and projectile combination, the incident energy and the angular momentum 
involved. In nucleon induced reactions, for example, the compound process is dominant in 
low energy region, while the pre-equilibrium and spallation processes become more likely 
as the incident energy increases. In heavy-ion collisions, we also have to introduce various 
models of different natures depending on each specific process. However, most of them are 
restricted to the specific energy regime or specific phenomenon and some of them have too 
many parameters to obtain a definite physical conclusion from the analysis. 

The main purpose of the series of our work is to develop a model which can describe the 
various aspects of nuclear reactions in an unified way. We try to seek a model with minimum 
number of parameters, a wide range of applicability, and a quantitative agreement with as 
many observables as possible. In addition to these requests, we demand the model to be so 
simple that one can run its computer code on work stations. 

In the heavy-ion physics, microscopic models, which describe reactions in terms of the 
dynamics of the interacting nucleons, are commonly used to extract the information of 
the nuclear matter under extreme conditions from the final observables. The most pop- 
ular models of this type are the Boltzmann-Uehling-Uhlenbeck/Vlasov-Uehling-Uhlenbeck 
(BUU/VUU) |], the quantum molecular dynamics (QMD) @, and CASCADE type mod- 
els Hd]. So far these microscopic models have shedded light on several exciting topics in 
heavy-ion physics, e.g. the multifragmentation, the flow of the nuclear matter, and the 
energetic particle production. However the parameters of the models such as the effective 
interaction, elastic and inelastic channels of NN cross section, differ substantially from one 
model to the other even in the same type of model. Furthermore, these models have not 
been tested intensively in much simpler light-ion reactions except for an analysis of (p,xn) 
reaction carried out by Peilert et al . In their analysis, however, the lower part of neutron 
energy spectra cannot be treated, since statistical decay following the QMD process was not 
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included. We thus start the series of our work from the analysis of the simplest type of the 
reactions, the (N, xN') (nucleon in, nucleon out) reaction in this paper, aiming to establish 
an unified model for various nuclear reactions. In the subsequent works, we are planning to 
analyze (N,xtt), (tt,xN) reactions and heavy-ion reactions. 

We restrict our subject to the reactions of nucleon-nucleus, meson-nucleus and nucleus- 
nucleus collisions with energies well above Coulomb barrier up to several GeV/nucleon, 
where the classical treatment of the collisions is justified in a first-order approximation. We 
do not deal with the phenomena which are dominated by the quantum effects. In this energy 
regime, the whole reaction process could be divided into two parts, i.e. dynamical process 
and statistical process. These two processes are well separated in their time scales. In the 
dynamical process, the direct reactions, nonequilibrium reactions, and dynamical formation 
of highly-excited fragments take place during typical collision times of the order 10~ 22 sec. 
After that, the evaporation and fission decay, which we call the statistical process, occur 
in the longer time scale of the order 10~ 21 — 10~ 15 sec. We thus employ a two step model, 
namely, we incorporate quantum molecular dynamics (QMD) for the dynamical process with 
statistical decay model (SDM). Similar hybrid models have been used in the analysis of the 
heavy- ion collisions 0-0]. 

In this paper, we define the basic ingredients of QMD plus SDM model and discuss how 
these two are combined without introducing any additional parameter. We then apply this 
model systematically to (N, xN') reactions, and discuss which element in the model is crucial 
for describing these reactions and what is necessary to develop the model further. In Sec. 
II we describe the details of the QMD, the effective interaction, NN elastic and inelastic 
cross sections, the relativistic corrections and the statistical decay model employed in our 
model. In Sec. Ill we compare the various double differential cross sections calculated by 
this model with the experimental data for proton induced reactions with incident energies 
from 100 MeV up to 3 GeV. We summarize and conclude this work in Sec. IV. 
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II. DESCRIPTION OF THE BASIC MODEL 



As we mentioned above, our basic model consists of two part, the quantum molecular 
dynamics (QMD) and the statistical decay model (SDM). The reason of employing the 
QMD model for the description of the dynamical processes is that the QMD can calculate 
the fragment formation in a natural and practical way Though the QMD method is widely 
used in the study of nuclear fragmentation 0, the details of the prescription differ from 
author to author. Aiming to establish a simple standard model, we will start from the 
standard type of QMD, taking into account of the relativistic kinematics and the relativistic 
correction for the effective interaction. Additionally, we treat the resonances of nucleon, 
A and A* (1440), and real pions with their isospin degrees of freedom in the equation of 
motion. For the statistical decay process, we use a simple prescription including only the 
light particle evaporation. 



A. Quantum Molecular Dynamics 



1. Basic equations and effective interaction 



The QMD method is a semi-classical simulation method in which each nucleon state is 
represented by a Gaussian wave function of width L, 
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where and Pj are the centers of position and momentum of i-th nucleon, respectively. 
The total wave function is assumed to be a direct product of these wave functions. Thus the 
one-body distribution function is obtained by the Wigner transform of the wave function, 
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The time evolution of R» and Pj is described by Newtonian equations and the stochastic 
two-body collision term. The Newtonian equations are derived on the basis of the time- 
dependent variational principle |2j as 

OH 



(4) 



dPi ' " ' <9R, L 

where the Hamiltonian H consists of the single-particle energy including the mass term and 
the energy of the two-body interaction. As for the effective interaction, we adopt the Skyrme 
type, the Coulomb, and the symmetry terms in this paper. By using the Gaussian function 
of nucleons [Eq. (|])], we get 
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with 



Ei = J ml + P, 2 



(6) 



where erf denotes the error function. In the above equation, q is 1 for protons and for 
neutrons, while < pi > is an overlap of density with other nucleons defined as 



with 



< Pi > =J2 P^j = H j dr Pi( r ) Pj( r ) 

= J2 (47rL)- 3/2 exp [-(R; - R^/AL 



(7) 
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In this paper we use the parameters A = —219.4 MeV, B = 165.3 MeV, and r = 4/3 which 
yield a compressibility of K = 237.3 MeV, saturation at p = p = 0.168fm -3 and a binding 
energy of 16 MeV per nucleon for infinite nuclear matter. The symmetry energy parameter 
C s is chosen to be 25 MeV. The width of Gaussian L is a parameter of the QMD and fixed 
as L = 2.0 fm 2 in this paper. 



2. Two-body collision term 

In addition to the Newtonian equation Eq. (f|), the time evolution of the system is 
affected by the two-body collision term. In the QMD method, the stochastic two-body 
collision process is introduced in a phenomenological way on the analogy of the test-particle 
calculation of the BUU collision term |l|. It includes the Pauli blocking factor (1 — /(r, p, £)), 
which is lacking in the cascade collision process 0,|J. We follow basically the prescription 
of the two-body collision term used in the BUU calculation done by Wolf et al. fl9,lQl, and 



modify it to extend the energy range up to 3 GeV. We thus describe here the outline of the 
procedure of Ref. PJTOfl and explain the extensions introduced in this paper. Further details 



of the collision term and the dynamics of A's, iV*'s, and pions discussed below can be found 
in Ref. 1,0. 

It is assumed that two particles collide if their impact parameter defined in a covariant 



way is smaller than a given value 6 max = yc/n obtained from cross section a. The collisions 
are considered as instantaneous interaction and a collision event is specified by the two points 
in space-time where the collision happens. Therefore it is hard to retain the covariance, since 
one has to choose a common reference frame for the QMD calculations. Hence the average 
proper time of the collision points defined by each particle is used to determine the time 
step in which the collision happens. This collision prescription was checked for heavy-ion 
collisions from 400 MeV/nucleon to 2.1 GeV/nucleon, and it was found that the disturbance 
of the covariance was very small ||. 

In order to treat the reactions with high bombarding energies, we include in our QMD 



6 



simulation the nucleons (N), deltas (A(1232)), AT*(1440)'s and pions with their isospin 
degree of freedom. The A's and iV*'s are propagated in the same interactions as the nucleons 
except for the symmetry term, while pions feel only the Coulomb interaction. The creation 
and absorption of these particles are treated in the collision term. In the following, we list all 
channels included in the collision term, where B denotes a baryon and N, more specifically, 
a nucleon: 
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(9) 



The channel 8 has been added to the prescription of Wolf et al. 0, which is the inverse 
process of the additional decay channel of A^*(1440)'s (cf. channel 11 below). 

We use the following parametrization for all baryon-baryon elastic cross section (channel 
1 in Eq. (I)), 

o = Cl _ +C 2 (mb), (10) 

1 + lOOv/i' V ; ' V ; 

with 

' = max(0, y/s-Mi-Mj- cutoff) (GeV), (11) 

where cutoff is 0.02 GeV for nucleon-nucleon channel, while it is zero for the others. This 
is the conventional Cugnon parametrization form We use this form up to y/s ' = 

0.4286 (GeV), which corresponds to 1 GeV Lab energy for nucleon-nucleon case. Above 1 
GeV, we parametrize the experimental p-p and p-n elastic cross section P~3"|.|T4" 



as, 



o = C, 



2 

l--tan _1 (l.5Vs'-0.8) +C 4 (mb). (12) 



In order to connect Eqs. ([H]) and (|1~2| ) smoothly, we slightly modified the parameters of 
Cugnon The actual values of the parameters C{ in the above equations are listed in 

Table |. The angular distribution of the elastic channels is taken from the same form as 
Cugnon parametrization JI|^|. 

For the production of baryonic resonances (channel 2 and 4 in Eq. (0)), we adopt the 
total cross section based on the parametrization of VerWest and Arndt fil l, in which the 



pion cross sections are parametrized assuming the pions are produced through baryonic 
resonances. Their parametrization was performed by fitting the experimental data up to 1.5 
GeV incident energy. In order to extend the energy range up to 3 GeV, we have modified 
the parameters in the following way. In the model of VerWest and Arndt, the cross sections 
are parametrized according to the initial and final total isospin % and / of the two-nucleon 
system |L1] as 

Tr(ftc) 2 (p r y mlY\ q /q Q f 
if{) 2tf [poj ( S *-m 2 ) 2 + m 2 r 2 ' K } 



where 
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s* = (M) 2 , g = ?K), 



and (M) is the mean mass of the resonance UTTJ obtained from a Breit-Wigner distribution 



with M = 1220 MeV, T = 120 MeV for the A and M = 1430 MeV and T = 200 MeV for 
the N*. From these cross sections, we determine the production cross section of A's [Of as 
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p + p -> n + A ++ : cr 10 + ±er u , 



p + p -> p + A+ : |o-ii, 



n + p -> p + A : \<j X \ + ±cr 10 , 

(15) 

n + p -> n + A+ : |crn + |<7i , 
n + n — > p + A" : a i + ^a n , 



n + n — > n + A : |o"n. 



We have effectively included the cross section of the nd final state, parametrized as af in 



IIJ , in the cross section of a w . 

We assume in this paper that the cross section cr 01 in [|ll|] contributes only to the N* 
production independently of the isospin components. Thus we rename aoi as &n*, and the 
iV* production cross sections are given by 
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The new parameters in Eq. ( |13D are given in Table |T[ In order to determine these parameters 
and the parameters of elastic cross section in high energy part defined in Eq. (fL2j), we fitted 
the experimental p-p and p-n cross sections |13|, 14 1 . In Fig. H we show the p-p (a) and p-n (b) 



total (solid line), elastic (long dashed line), and inelastic (dot-dashed line) cross sections. 
The inelastic cross section is the sum of the A (short dashed line) and iV* (dotted line) 
production cross section, calculated by Eqs. (|1^,|1^,|1^,|1^,|1^) . In the same figure, we show 
the corresponding experimental total (open circles), elastic (open triangles), and inelastic 
(open boxes) cross sections JTJJ with error bars. For the p-p case, the present parametrization 



of elastic and inelastic cross sections fits to the data for whole energy range up to 3 GeV 
except for the some deviation around 1 GeV, which is due to the elastic cross section and 
does not affect the result. On the other hand, for p-n case where only the total cross section 



is available in the data, we fitted it at the energy higher than 0.7 GeV up to 3 GeV, and 
adopt the Cugnon's type elastic cross section in the low momentum region instead of the 
free elastic cross section. 

In Fig. ^ we show the pion cross section of pn — > nnir + + pp7t~ (solid line) obtained by 
our new parametrization of Table [TT]. In the same figure, the gray bold line denotes the result 
of the original parametrization of VerWest and Arndt [pi , while the experimental data are 



taken from Ref. |14fl . By this new parameter set, our pion production cross section below 1.5 



GeV does not differ so much from the original results of VerWest and Arndt (TT| , which are 



essentially the same as the data. However, above 1.5 GeV, our result fitted the experimental 
data, while the result obtained by extrapolating the original parametrization of VerWest 
and Arndt to higher energy shows a big bump, which has no experimental support. 

In the higher energy region, the role of N* becomes important. One of the good quan- 
tities which shows the characteristics of the higher resonances is the elementary two pion 
production cross section. In the present prescription it is described only by a^* combined 
with the decay modes of the resonances which will be mentioned below. For example, it is 
shown in our prescription that 

a(pn -> ppir n-) = -j-(T N *, 
15 

5 

a(pn -> pnTt + 7r~) = ^n*, (17) 
a(pp -> ppix + ix~) = -a N *. 

We thus plot the a^* in Fig. |3| as well as the experimental two pion production cross 
sections Jl4| scaled by above factors. This figure shows that our parametrization fitted the 
gross features of the experimental data for the energy range up to 3 GeV. Although this 
parametrization should be modified if two delta production channel or the direct two pion 
decay of N* or higher resonances are included, the present prescription of the elementary 
inelastic channels could roughly reproduce the experimental single and two pion production 
cross sections for the energy range up to 3 GeV. 

We do not take into account the direct s-state pion production mechanism but all pions 
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are assumed to be produced through baryonic resonances. The masses of the resonances 
are randomly distributed according to the Breit-Wigner distribution with the momentum- 
dependent width IJI^I , i.e., 

25r 2 

f M = (M-M r )* + 0.2^ ' < 18) 



with 

' q \ 3 M r (v( q )Y T ^ (ig) 



\q r J M \v(q r 

where q denotes the cm. momentum in the irN channel, and index r refers to the values at 
the mass M r , and 

= ^ (20) 

We have applied this momentum dependent width not only to the A-resonance but also to 
iV*(1440). The values of the parameters used in this paper are listed in Table |T|. 

Another important ingredient of the resonance production (channel 2 and 4 in Eq. (^)) 
is the angular distribution of the resonances. Wolf et al. || parametrized the angular 



distribution of the experimental data ]T9| for p + p — > n + A ++ and assumed the same 



angular dependence for each isospin channel in the following form, 

<7r(s, cos#) — a (s) {ai(s) + 3a 3 (s) cos 2 9], (21) 

with 

ao(s) = , , , \, , vr - (22) 
47r(ai(s) + a 3 (s)) 



The values of a\(s) and a^(s) are given in Table [IV|. In the high energy region y/s > 2.4 GeV, 
which corresponds to the laboratory energy higher than 1.2 GeV, this angular distribution 
is assumed to be constant, since there is no experimental data to be fitted in this energy 
region. However, above the resonance region E^b > 1-2 GeV, this is not justified because, for 
example, the angular distribution of protons from 27 Al(p,p') at 3.17 GeV reaction calculated 
by Eq. (|2l|) deviates from the experimental data (cf. discussions on Fig. § in the next section). 
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In order to get a better parametrization for high energy part, we assume that the angular 
dependence is effectively written as a sum of g R and another term g D as 



where 



g A (s,cos9) = - [g R {s, cos 9) + g D {s, cos 9)}. 



(23) 
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(27) 



This form of g^ is obtained by modifying the Cugnon parametrization |T|Q] of the NN elastic 
angular distribution so as to trace the angular distribution of Eq. ([H]) in the resonance region, 
and approach the elastic-like angular distribution for the higher momentum region. We use 
this angular distribution for A-resonance and iV*(1440). For the latter case, we apply this 
formula by shifting the energy ^/s by the mass difference of the two resonances, i.e., 208 
MeV. The energy dependence of this angular distribution is shown in Fig. ^, where we plot 
#R (gray bold dashed lines) and <?d (solid lines) in (a), while #r (gray bold dashed lines) and 
gx (solid lines) in (b). In these figures, we symmetrized the elastic-like angular distribution, 
i.e., ![#d( s j cos9) + gr>(s, — cos9)}, in order to compare it with the angular distribution of 
Eq. fl2l| ) (gray bold dashed lines). In the next section, we will discuss the dependence of the 
final results of (N, xN') reactions on this angular distribution of the hadronic resonances. 

The cross sections for the channel 3 and 5 in Eq. (|^) are determined by the law of 
detailed balance from the cross section of channel 2 and 4 with taking into account the mass 



dependence of the cross section fL0|| . 
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For the pion absorption cross section on nucleons (channel 6 and 7 in Eq. @) we take 
the maximum cross sections from the particle data table |TJ| and scale them according to 
the Breit-Wigner formula. For the case of pion absorption on A (channel 8 in Eq. (|9])), 
we assume the same cross section as the channel 7 with shifting the energy by the mass 
difference of the A and nucleon. 

Apart from the collision term, we take into account the decay of the baryonic resonances 
during the propagation as 

9. A ->• iV + vr 

10. N* -> iV + 7T (28) 

11. N* — > A + 7T 

The decay probability of the resonances is determined by a exponential decay law using 
their momentum-dependent width Eq. ([19]) and their proper time. The decay is assumed 
to be isotropic in their rest frame. The branching ratio of the channel 10. and 11. of Eq. 



is taken from the particle data table |L3| as T n *_>a +7T /[T n *_>a +7T + rjv.^jv +7r ] = 0.4. 



The other branching ratios concerning their isospin are determined from the appropriate 
Clebsch-Gordan coefficients. 



3. Relativistic corrections 

Non-covariant framework is another problem of applying the QMD method to the reac- 
tions at higher bombarding energies. As explained above, we have already introduced the 
relativistic form of energy expression in the Hamiltonian Eq. (j5|) and adopted the relativistic 
kinematics in the collision term. However a covariant formulation of the interaction term 
is necessary for a full relativistic description. A Lorentz covariant extension of the QMD, 
dubbed relativistic quantum molecular dynamics (RQMD), has been proposed by Sorge et 
al. |T5| based on the Poincare-invariant constrained Hamiltonian dynamics. Although the 
RQMD is a numerically feasible extension of QMD toward a fully covariant approach, it 
still costs too much computing time to apply the RQMD model to heavy systems. We 
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thus make the following alternative extension of QMD and include the main part of the 
relativistic dynamical effects in our model. 

Lehmann et al. |IB] have compared the time evolution of the phase space and the particle 
production obtained by QMD and RQMD, looking for the relativistic effects in heavy-ion 
collisions in the intermediate energy regime. They found that there is no significant difference 
between the results of QMD and RQMD in the r] and 7r meson production cross sections 
j~6| and the proton inclusive spectra . The difference appeared only in the values of the 



maximum density |l7j and the transverse flow Both are larger in RQMD. Their studies 



showed that a large part of the difference comes from the Lorentz contraction of the initial 
phase space distribution in RQMD. If this Lorentz contraction is employed in the normal 
QMD, however, the transverse flow is overestimated. Thus they concluded that in RQMD 
this effect is partially counterbalanced by the covariant treatment of the interaction, but 
there still remains an increased flow compared with the normal QMD calculation. Based 
on their investigation, we introduce in this paper the Lorentz-scalar quantities into the 
arguments of the interactions in Eq. (^j) as well as the Lorentz contraction of the initial 
phase space distribution. By these modifications the main part of the relativistic dynamical 
effects would be approximately included in our QMD. 

All arguments of the interaction in Eq. (y) are written as a function of the squared spatial 
distance: 

R| = (R, - R,) 2 . (29) 



In the RQMD ||16|jl7| , these arguments are replaced by the squared transverse four- 



dimensional distance — q% .. defined as 

= -4 + (30) 

Pij 



where qij is the four dimensional distance g« — Qj, while pij denotes the sum of the four 
momenta of the two particle pi + pj. In the c.m.s of the particle % and j, the squared 
covariant transverse distance —q^. reduces to the usual squared distance. We therefore 
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change the argument in Eq. (^) from the Rf • to the squared distance in the c.m.s. of the two 
particles R?-, where the tilde means the quantities defined in the c.m.s of the two particles, 

RJ = HJ+7S(IV&) 2 , (31) 

with 

^=tt^- (32) 

By this change, all interactions of the Hamiltonian [Eq. (|5|)] depend also on the momen- 
tum. The form of the equation of motion [Eq. (|4])] changes to 

R^tI^ + EA^, (33) 



m 2 + p2 Y dP 
> 

9R 



(34) 



.i 



with 



IA 1 

Dij ~~27 2L Pij 



1 BT, T _! , r _ x N 1 



+ ^ C 4~E erf (^) +MP 4Rl 



I? I \ V ^/ , 



and 



C I 

(1 — 2|cf — c,-|) 2l Pij, (35) 



|5| = 2Ry + 2 7 | (R ?J • 4) 4 (36) 
E i + E j K^'^J 

x |r<* + il K • 4) (4 - §) } • ( 3? ) 

where is defined in Eq. (^). 

We also introduce the Lorentz-scalar quantities into the one-body phase-space distribu- 
tion function Eq. @ as 
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where Pf • denotes the squared relative momentum in the c.m.s of the particle i and j, which 
is expressed for two particles with the same mass as 

PJ = Pj - {E, - Ej) 2 , (39) 

with 

P« = P, - Pi- (40) 

At the starting point of the QMD calculation, we boost the ground state of projectile 
(and target as well if c.m.s of the target and projectile is chosen as a reference frame) 
according to the beam energy. The coordinate and momentum P]? z of the nucleon in 
the beam direction z after the boost are obtained by the Lorentz transformation from those 
before the boost, H iz and P iz , as 

R^ = (R J2 -R ,)/7 + Ro„ (41) 
PL = 7 (Pi* + (42) 

where Ro denotes the initial cm. coordinate of the nucleus, while (3 and 7 are the boost- 
ing velocity and its gamma factor, respectively. At this moment, the potential energy of 
the system and the phase-space distribution function keep the same values as those before 
the boost due to their Lorentz-scalar properties discussed above. During the propagation 
of the boosted nuclei, however, those quantities are changing slightly even in the above 
prescriptions, since the equations of motion Eq. ([33|) and (|34|) are not covariant. But the 
disturbance due to the non-covariant feature of the equation of motion is negligibly small up 
to the energy 3 GeV/nucleon. We thus introduce the relativistic corrections discussed above 
to the non-covariant QMD to save the computing time instead of using the full covariant 
framework. 

We should mention here that if one employs the Lorentz contraction for the boosted initial 
state but does not replace the arguments of the interaction and phase-space distribution by 
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the Lorentz-scalar ones, the potential energy decreases about 80 MeV and the phase-space 
factor at each nucleon's point changes about 40 % after the boost in the case of 40 Ca even at 
1 GeV/nucleon boosting energy. This means that the boosted ground state obtained by this 
way may decay spontaneously before it collides with the other nucleus. By our prescription, 
we are free from this problem. 

We have checked the above prescription by analyzing the transverse flow, which is sensi- 
tive to the treatment of the relativity as discussed in Ref. |R|. In Fig. || we show the energy 



dependence of so called directed transverse momentum < P^ ir >, which is a measure for the 
transverse flow and defined by 

1 N 

< p f >= ITr E si g n [ y CM«] (43) 
iV i=i 

where 1cm {i) is the rapidity of the z-th baryon in the c.m.s and its transverse momentum 
in the reaction plane. We plot the result of the present QMD simulation (solid line with 
full boxes) as a differences from that of RQMD [IjJ for 40 Ca + 40 Ca reactions at b = 2 fm, 
for the energy range from 150 MeV/nucleon to 4 GeV/nucleon. In these calculations, we 
use the same ground states (mentioned below), the same gaussian wave packets and the 
same interaction (mentioned before). We omitted the Coulomb interaction and two-body 
collision term for a simplicity. In Fig. [|, we also plot the other results obtained by the 
standard QMD (dot-dashed line with open circles) without the initial Lorentz contraction 
and without the relativistic corrections, and by the standard QMD only with the Lorentz 
contraction (dashed line with open boxes). As mentioned before, the Lorentz contraction of 
the initial phase space distribution increases the flow, which is shown by the change from 
the dot-dashed line to the dashed line. By the full covariant treatment of the interaction, 



however, this effect is counterbalanced, but still remains an increased flow |l6fl . As seen in 
Fig. [5], our prescription does not deviate so much from the full covariant treatment up to 3 
GeV/nucleon. We thus expect that our QMD simulation with the relativistic corrections is 
very close to the covariant simulation RQMD in this energy regime. 

At much higher energy, however, our result is decreasing linearly from that of RQMD. 
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This deviation comes from the different treatment of the potential between the RQMD fl7 
and our QMD; a Lorentz scalar type in the former, while a time-component of the vector type 
in the later, respectively. This is understood qualitatively by considering a single particle 
motion under a fixed external potential U. In the Lorentz scalar treatment of the potential 
U, the single particle energy p® is expressed in this simple case as 



Pi = yPi + m i + 2m i u - ( 44 ) 
Accordingly the equation of motion is 

P. = (45) 
Pi dRi 

On the other hand in our prescription they are 



Pi = \jPi + mf + U, (46) 

and 

- 9U (47) 



In this test calculations, the form of J|r-, which is attractive in the beginning, is almost 
the same in our QMD and in the RQMD. Thus the force of our QMD becomes larger 
and deviates linearly from that of the RQMD as energy increases. Above 3 GeV/nucleon, 
therefore, the full covariant prescription is necessary to describe the reactions particularly 
the nucleus-nucleus collisions. This is out of scope of this paper. 

Some details of actual numerical calculations should be mentioned here, since all potential 
terms depend on the momentum by the relativistic corrections. In order to keep a numerical 
accuracy, we use the second order Runge-Kutta method to integrate the equations of motion. 
For the energy conservation for the collision term, we assume 

Ei + Ej + t/pot = E[ + E> + f/ pot , (48) 

where Ei, Ej, and E[, Ej are the energies of particle i and j before and after the collision, 
while U pot and U' t the potential energy of the system. We determine iteratively the final 



momenta of the colliding particles so as to satisfy the energy conservation Eq. (f48|). This 
prescription is applied to the channels 1 to 5 in Eq. (|9|), 9 to 11 in Eq. fl28|). For the pion 
absorption channels of 6, 7 and 8 in Eq. (|9|), the energy conservation is written as 

E i + E j + U pot = E R + Ul ot , (49) 

where Er is the resonance energy. In this case, we determine iteratively the rest mass of 
the resonance to conserve the energy. 

4- Properties of "ground state" 

An important ingredient of the QMD calculation is how to determine the initial phase 
space distribution of the projectile and target. For that we cannot use the real ground 
state (energy minimum state) of the system defined by the Hamiltonian Eq. (||p, since the 
model does not have the Fermionic properties. However, it is necessary to obtain a stable 
"ground state" . Also some typical properties of the real ground state should approximately 
be fulfilled, especially the binding energy and phase space distribution. To get such "ground 
state", we employ the following random packing procedure ||; 

We distribute the centers of position Rj of the individual Gaussian wave packet according 
to a distribution of the Woods-Saxon shape with the radius Rq = 1.124 A 1 / 3 — 0.5 fm and 
the diffuseness parameter a = 0.2 fm. We cut off the Woods-Saxon tail at -R max = 1.124 A 1 / 3 
fm. In this procedure, we impose minimum distance between the centers of the Gaussians 
in order to reduce the density fluctuation. We use 1.5 fm for the identical nucleons and 1.0 
fm for the other. 

Now we can calculate the density and potential energy at any point (here, we do not 
need the relativistic correction discussed in the previous sub-section). Then the center of 
momentum Pj is randomly sampled from the sphere of radius pp(Ri) which is the Fermi 
momentum obtained by the local Thomas-Fermi approximation. This sampling is rejected 
and another value is sampled if the sum of kinetic and potential energy of the particle is 
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positive or the phase space factor /(Rj,Pj) (cf. Eqs. (0) and (|3|) ) for the nucleon j which 
have been previously accepted violates the Pauli principle H. 



Finally, we check the total binding energy with the simple mass formula |20|| , i.e., 

E hin = - 15.56 A + 17.23 A 2/3 

+ 46.57 <^£ + ?-?£L (50) 
2 A 5 1.24 A 1 / 3 v 7 

If the binding energy per nucleon obtained by our sampling lies within Eu n ± 0.5, we adopt 
this configuration "ground state" . 

Thus the "ground state" obtained by this procedure always has an appropriate binding 
energy. However, there is still open phase space below the Fermi surface, since the "ground 
state" is not the energy minimum state of the Fermions. In fact, during the time evolution of 
the "ground state" under the QMD dynamics described in previous sub-section, only 70 % 
of the collisions are blocked by the final state Pauli blocking in the two-body collision term. 
It is allowed at a collision that one nucleon goes down to the lower energy state and the other 
goes up to the positive energy state. This means that some nucleons could be spontaneously 
emitted after some time due to the fluctuation of the configuration. To avoid this problem, 
we assume from a technical point of view that any pair of nucleons originated in the same 
nucleus does not collide each other until at least one of them experience a collision with a 
nucleon from the other nucleus. By this assumption the number of the emitted nucleons 
from the "ground state" is reduced to less than about 1 % of the nucleons up to the time 
200 fm/c. 

The density profile of the "ground state" obtained here has high density in the center 
and rather wide surface shape. This is due to the large width of Gaussian L = 2.0 fm 2 
used in this paper. On the other hand, the momentum distribution of the "ground state" 
well reproduces the result of Hartree-Fock calculation. In Fig. |6], we show (a) the density 
p(r) distribution and (b) the momentum distribution p(p) for the "ground state" of 40 Ca 
obtained by QMD simulation (solid lines). The results shown here are averaged quantities 
over time evolution up to 200 fm/c and over 100 events. The error bars in Fig. |6] denote 
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the fluctuations in time evolution averaged over 100 events. Although the fluctuation of the 
one event is much larger, this figure shows that the ground state profile is very stable in 
time on the average over the events. In the same figure, we also plot the results of Hartree- 
Fock calculation (dot-dashed lines) and the limit of infinite nuclear matter (dashed lines). 
The energy spectra of the emitted particles given in the next section, particularly of the 
subthreshold particle production [^TJ are more sensitive to the momentum distribution than 
the density profile. This is the reason why we adopted a parametrization which leads to a 
better momentum distribution at the cost of diffuse density profile. 



B. Statistical Decay Model 

At the end of the dynamical stage of the reactions, the QMD simulation yields many 
fragments, which are normally in highly excited states. One may think that the decay process 
of the excited fragments might be described by the QMD dynamics if we can continue the 
calculation for enough long time. However, we do not follow this method but instead we 
stop the QMD calculation and switch to the statistical decay model (SDM) at the end of 
the dynamical stage. There are two reasons for this hybrid model. One is that the time 
scales of the dynamical and statistical processes are quite different. It is not clever or even 
practically impossible to continue the reliable QMD calculation more than 10~ 20 sec which 
is necessary to calculate the decay process. Another is a more fundamental reason that the 
Fermi statistics, which is essential to describe the decay process of the fragments, cannot be 



traced correctly in the QMD simulation [22 



We identify the fragments together with their excitation energies at about 100 ~ 150 
fm/c of the QMD simulation. The dependence of the final results on this switching time 
will be discussed in the next section. Each fragment is recognized by using a minimum 
distance chain procedure, i.e., two nucleons are considered to be bound in a fragment if the 
distance between their centroids is smaller than 4 fm. We then calculate the total energy of 
the fragment in its rest frame and estimate the excitation energy by subtracting the ground 
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state energy given by Eq. (|50|). 

Though there have been proposed many sophisticated statistical decay codes so far, we 
use here the simple model of light particle evaporation. We consider only n, p, d, t, 3 He, 
and a evaporation. The emission probability P x of these particles x is given with the Fermi 
gas model as, 

P x = {2J x + l)m x ea x (e)p(E)de, (51) 

where J x ,m x , and e are the spin, mass and kinetic energy of the particle x, while a x (e) and 
p(E) denote the inverse cross section for the absorption of the particle with energy e and 
the level density of the residual nucleus with the excitation energy E, respectively. We use 
the following simple form for p(E), 

p(E) = w exp (2v / oE) , (52) 

with a = A/8 MeV" 1 and Wo is a constant. The inverse cross section is assumed to have the 
form 

[ (1 - U x /e) tcR 2 : e > U x 
a x {e) = I (53) 

jo : e < U x 

where R denotes the absorption radius and U x is a Coulomb barrier for the particle x, 



for which we employ empirical values used in the existing statistical decay code |23| . The 
excitation energy E in Eq. (|52|) is given by 

E = E -e-Q, (54) 

where E denotes the excitation energy of the parent nucleus and Q is the reaction Q-value 
calculated from the mass formula Eq. fl5"0|). The total emission probability R x of the particle 
x is obtained by integrating the available energy of Eq. (|5T]) as 

Rx = (2 J x + 1) m x 

rEo—Qx 

x ea x (e)p(E -Q x -e)de. (55) 

J U x 
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This integration can be calculated analytically and the energy spectrum of the emitted 
particles is given by 

N(e x ) de x = — ^ x expj- ^^^' j de x , (56) 

with 

aT* = E -U x -Q x . (57) 

In this formulation, we do not consider the 7 decay nor the angular momentum de- 
pendence. The latter is important for the heavy-ion reactions but not so serious for the 
nucleon-induced reactions considered in this paper. We simulate the whole statistical de- 
cay process as a sequential light particle evaporation discussed above by making use of the 
Monte-Carlo method until no more particle can be emitted. 



III. ANALYSIS OF THE (N,XN') REACTIONS 

In this section, we systematically apply the QMD plus SDM method described in the 
previous section to (N, xN') (nucleon in, nucleon out) type reactions. In order to get suffi- 
cient statistics, we performed the QMD calculations for a large number of events, typically 
50000 events, and averaged them to obtain the following results. 

We first check the dependence on the switching time t sw when the QMD calculation 
is stopped and switched to SDM, which is an ambiguous point of the present model. In 
Fig. |7|, we show a typical neutron energy spectrum at 30° laboratory angle for the reac- 
tion p(1.5GeV)+ 208 Pb. Note that x-axis is plotted in a logarithmic scale to compare the 
calculated results in detail with the experimental data particularly in a low energy region. 



The experimental data (full boxes with error bars) are taken from Ref. The solid his- 
togram denotes the final result of the QMD + SDM calculation. In this case, we switch the 
QMD calculation to SDM at 100 fm/c. In the same figure, we also plot the spectrum of 
the neutron obtained only by the QMD calculation up to 100 fm/c (dot-dashed histogram) 
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and that coming from the QMD fragments calculated with SDM (dashed histogram), re- 
spectively. The former shows a cascade and/or preequilibrium energy spectrum, while the 
latter an evaporation spectrum. These two components of spectrum are affected by chang- 
ing the switching time £ sw . However, the total spectrum shape, which is the sum of the two 
components, stays almost unchanged if an adequately long time is chosen for the switching 
time t sw . This is shown in Fig. H, where we plot results of the total spectra calculated by 
QMD + SDM with three different switching times, 50 fm/c (dashed line), 100 fm/c (solid 
line) and 150 fm/c (dot-dashed line). This figure shows that although the latter two lines 
resemble each other, they deviate definitely from the first line. This indicates that the QMD 
fragments before 100 fm/c are not in thermal equilibrium and that within a time interval 
from 100 fm/c to 150 fm/c the decay processes of the excited fragments described by QMD 
and SDM are nearly equivalent. Although we should keep in mind that these two are not 
identical at low temperature as the former is always dominated by the classical statistics 



22| , we can conclude that the final results are not sensitive to the switching time £ sw as long 
as it is chosen after the time when the thermal equilibrium is achieved and before the time 
the temperature of the fragments become low and classical statistics breaks down seriously. 
The similar conclusion has been obtained in Ref. ||, which indicates that the minimum 
switching time to get stable results depends on the size of system and the incident energy. 
For the case of nucleon induced reactions, we found that 100 fm/c is enough to get stable 
results and we use this value for all systems in the present study. 

Next check is the detailed examination of the inelastic channels in the two-body collision 
term. For this purpose we compare our results with the experimental data at high incident 
energy of proton on the light-mass target, which directly reflects the elementary processes 
included in the model. In Fig. |9|, we plot the invariant cross section of the proton (left-hand- 
side) and negative pion (right-hand-side) emission for the reaction p (3. 17GeV)+ 27 Al. The 



experimental data (full boxes with error bars) are taken from Ref. [^] and the results of 
QMD + SDM are denoted by solid histograms. In the same figure, we plot results of QMD 
+ SDM with the different choice of the angular distribution of the resonances. The dashed 
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histograms are the results obtained with only the resonance-like angular distribution of Eq. 
(pl|), while the dot-dashed histograms are those with the direct-like angular distribution of 
Eq. fl2"l|), respectively. This figure shows that the average of the two components of the 
angular distribution of Eq. ( |23|) well fits the experimental proton spectra. On the other 
hand, the pion spectra are rather insensitive to the angular distribution of the resonances. 
Instead, their spectra are dominantly determined by the mass distribution of the resonances 
of Eq. (|18|). Although the authors of Ref. p5j analyzed these data by making use of the 



two-moving-source model, our QMD + SDM can reproduce excellently the proton and pion 
spectra at the same time without any special assumption. 



In Figs. |10|-|T5|, we compare the neutron energy spectra obtained by the QMD + SDM 
calculations with the experimental data p4] , |2~6|| for Fe and Pb targets at the proton energy 
from 113 MeV up to 3 GeV. In the fields of application of accelerators, such as spallation 
neutron sources, accelerator-based transmutation systems, and shielding of cosmic ray in 
space activity, the production of slow neutrons plays an important role. That is the reason 
we chose these data [^3|,^BJ to compare with, since the neutron spectra from 1 MeV up to 
the beam energy are available. For efficient comparisons of the calculations and the data, 
both for low energy and high energy regions, we plot the same results in two figures with 
the x-axis in a logarithmic scale (left-hand-side) and in a linear scale (right-hand-side). One 
can see the detail of the thermal and preequilibrium neutron spectra in the former figure, 
while the direct-like components of the spectra in the latter figure, respectively. 

These figures indicate that over the broad range of the incident energies from 100 MeV to 
3 GeV, independently of the targets, and of all angles of the outgoing neutrons, our results 
of the neutron energy spectra agree well with the data from 1 MeV up to the beam energy. 
Though one may notice some disagreement at the high energy part of the most forward 
angle, which will be discussed later, the overall agreement is satisfactory. Particularly a 
remarkable agreement of the present calculations with the low energy neutron data below 
several tens MeV confirms that the QMD gives proper excitation spectra of the excited 
fragments from which the statistical neutron emission takes place. With a suitable chosen 

25 



fixed set of parameters, our QMD plus SDM model is able to reproduce quantitatively the 
overall neutron spectra for the broad range of the incident energies and target masses. 

At 113 MeV, we additionally compare our results with the prediction of the intranuclear 
cascade plus light particle evaporation model (NUCLEUS) P?fl . This model is essentially 
the intranuclear reaction part of NMTC |2S| and HETC codes. Calculations with 



NUCLEUS yield almost the same results as ours in the forward angles but give lower values 
in the backward angle, which is denoted by the dashed histogram at 150 degree in Fig. pi]. In 
this energy regime it has been reported [[3(],[nj that the semiclassical preequilibrium models 
based on an intranuclear cascade model also fail to reproduce the angular distributions. 
We found from the detailed comparison of the calculations that the underestimation of the 
backward angle in the above models is due to the insufficient treatment of the soft interaction 
of a nucleon with all the rest of the nucleons in the nucleus, which is naturally included in 
the QMD formalism. 

One may think that this explanation is in contradiction with Ref. PJ33|| where they at- 



tribute the failure to the insufficient contributions from second and higher order collisions. 
To resolve this problem, we have checked that the NUCLEUS has almost the same prescrip- 
tion of hard nucleon-nucleon interaction and has almost the same momentum distribution in 
the ground state as QMD. Difference is that our QMD part includes the soft nucleon-nucleon 
interaction but NUCLEUS does not. This soft interaction diffracts the nucleon. As a result, 
the yields of the backward angle increase and by the same reason the number of multiple 
hard collisions also increases. 

On the other hand, a multistep model of the Feshbach-Kerman-Koonin (FKK) |32| has 



been also applied to the energy regime around 100 MeV [p^] . Although the FKK approach 
successively reproduced the angular distribution, the overall absolute values of their results 
are very sensitive to the strength parameter of the residual interaction, which is adjusted 
to fit the experimental data. The strength so determined depends on the incident particle, 
target nucleus and incident energy. In QMD, on the contrary, the parameters of the soft 
nucleon-nucleon interaction in Eq. (g) are taken common to all reactions and determined 
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from the nuclear saturation condition. In addition, the final results are not so sensitive to 
them. 

The first analysis of (p, xn) reactions by the QMD approach in the energy regime up to 
800 MeV has been done by Peilert et al. ||. The neutron spectra of their results are very 
similar to those of the present work above several tens MeV. Their analysis, however, cannot 
predict a whole spectra of neutron, since the contribution of the statistical decay from the 
excited fragments produced in the QMD calculation was not considered in their work. 

Though the present results show overall agreement with the data for the very broad 
energy regime, one can see a systematic deviation from the data in the high energy part 
of the neutron spectra at the most forward angle at incident energy from 113 MeV up to 



800 MeV (see the right-hand-side of Fig. |T0| - |r3"|) . We suppose that the soft nucleon-nucleon 
interaction is responsible for this deviation. One possibility is a momentum dependent 
interaction that is not included in the present QMD, by which the nucleon could be affected 
coherently by the surrounding nucleons when its momentum is drastically changed by the 
hard nucleon-nucleon scattering. For the higher incident energies (see Fig. |TJ] and |T3p, this 
deviation disappears. In those cases, we have checked that the neutrons in the high energy 
part of the forward angle emerge after at least once experiencing the resonances of nucleon, 
and that the effects of the soft interaction is relatively small. The analysis by the QMD 
including the momentum dependent interaction will be reported in the forthcoming paper. 



IV. SUMMARY AND CONCLUSION 

We have proposed the quantum molecular dynamics (QMD) incorporated with the sta- 
tistical decay model (SDM) aiming to describe various nuclear reactions in a unified way, 
and applied this model to the (N, xN') reactions. We have checked and found that the 
final results do not depend on the switching time when the QMD simulation is stopped and 
switched to the SDM calculation as long as the switching time is chosen between 100 fm/c 
and 150 fm/c for the nucleon induced reactions. Therefore there left little ambiguity with 
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respect to the switching of two different kinds of models to describe whole process of the 
reactions in a unified way. 

In order to describe the reactions at high incident energies up to 3 GeV, we have taken 
into account two baryonic resonances, the A(1232) and the iV*(1440) as well as the pions in 
the QMD model. The elementary cross sections related to these resonances and pions are 
basically taken from the experimental data. The angular distributions of the resonances, 
which information is very poor in the experimental data, have been fixed to fit the 27 A\(p,p') 
data |25j at 3.17 GeV. It should be noted that the energy spectra of the nucleons from 
the (N, xN') reactions on the small target are suitable quantities to obtain the detailed 
information of the angular distribution of the resonances, while the pion spectra are useful 
to extract the information of the mass distribution of the resonances. 

In addition to the relativistic kinematics and approximately covariant prescription of 
the collision term, we have introduced the Lorentz-scalar quantities to the arguments of the 
interactions and to the phase space factor. By this relativistic corrections together with the 
Lorentz contraction of the initial phase space distribution, the main part of the relativistic 
dynamical effects is approximately described in our QMD for the energy regime up to 3 
GeV/nucleon. Validity of this model has been confirmed by the analysis of the transverse 
flow for the heavy-ion collisions in comparison with the results obtained by the covariant 
version of quantum molecular dynamics (RQMD). 

We have applied systematically QMD + SDM to the (N, xN') reactions for a broad 
range of incident energies from 100 MeV to 3 GeV and of target masses. The present 
model reproduced the overall feature of the outgoing neutrons quite well without assuming 
any reaction mechanism, and without changing a parameter set. Although there are a lot of 
parameters in the model which have not been investigated extensively in this paper, the final 
neutron spectra analyzed here do not depend so much on them, for example, the equation 
of state (choice of the interaction), width of the gaussian wave packets, and the detail of 
the statistical decay process. The main ingredients of the model, which produce the present 
results of the neutron spectra down to the energy of several MeV, are the parametrization of 
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the elastic and inelastic elementary cross section and the many-body dynamics itself, which 
have been discussed both in detail in this paper. We thus conclude that the present QMD 
+ SDM scheme gives a unified picture of the major three reaction mechanisms of (N,xN') 
reactions; i.e. compound, pre-equilibrium and spallation processes. 

Finally, we should mention that the present model is ready to be applied directly to the 
heavy-ion reactions in its original form. A study of heavy-ion reaction using this model is 
now under consideration. 
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FIGURES 

FIG. 1. p-p (a) and p-n (b) total (solid line), elastic (long dashed line), and inelastic (dot-dashed 
line) cross sections, which is the sum of the A (short dashed line) and N* (dotted line) production 
cross section, calculated by Eqs. (|i~0| , |l2[ |l~3| , [l5| , [l~6|) . The experimental total (open circles), elastic 
(open triangles) and inelastic (open boxes) cross sections are taken from Ref. Jl4j. 

FIG. 2. Pion cross sections of pn — > nnir + + ppir~ (solid line). The gray bold lines denote the 
results of the original parametrization of VerWest and Arndt pi , while the experimental data are 



taken from Ref. 14]. 



FIG. 3. Two pion production cross section. The solid line denotes ctn*, while the experimental 



data [14 1 are scaled according to the factors of Eq. (|i7|). 



FIG. 4. The energy dependent angular distribution of resonance production, (a) Two compo- 
nents of the angular distribution. The solid lines represent the elastic like angular distributions 
of Eq. (|24]), |[go(s, cos 9) + gr>{s, — cos 9)], while the gray bold dashed lines are the resonance like 
angular distributions of Eq. (pl|). (b) The total angular distribution of A (solid lines), which is the 
sum of 5r and go- The gray bold dashed lines are the same as in (a) 

FIG. 5. The directed transverse momentum as a function of energy /nucleon for 40 Ca + 40 Ca 
reaction at b = 2 fm. The results are shown as the differences from that of the RQMD [|l^| . The 
solid line with full boxes denotes the result of the present model with the relativistic corrections 
and the initial Lorentz contraction. The dot-dashed line with open circles is the result of the 
standard QMD without the initial Lorentz contraction and without the relativistic corrections, 
while the dashed line with open boxes is the result of the standard QMD only with the initial 
Lorentz contraction, respectively. 
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FIG. 6. (a) Density distribution p{r) and (b) momentum distribution p(p) of the ground state 
of 40 Ca obtained by the QMD simulation (solid lines) averaged over the time evolution up to 200 
fm/c and over 100 events. The error bars denote the fluctuations in time evolution averaged over 
100 events. The dot-dashed lines show the results of Hartree-Fock calculations, while the dashed 
lines denote the limit of infinite nuclear matter. 

FIG. 7. Neutron energy spectrum at the 30° laboratory angle for the reaction 
p(1.5GeV)+ 208 Pb. The experimental data (full boxes with error bars) are taken from Ref. [24]. 



The solid line denotes the final result of the QMD + SDM calculation with the switching time 
i sw = 100 fm/c. The dot-dashed line is the result obtained only by the QMD calculation up to 100 
fm/c, while the dashed line the neutron spectra coming from the QMD fragments calculated with 
SDM, respectively. 

FIG. 8. The total energy spectra calculated by QMD + SDM with different switching time, 
tsw = 50 fm/c (dashed line), i sw = 100 fm/c (solid line) and t sw = 150 fm/c (dot-dashed line). The 
reaction system is the same as in Fig. [?]. 

FIG. 9. Invariant cross sections of the proton (left-hand-side) and negative pion 
(light-hand-side) emission for the reaction p(3.17GeV)+ 27 Al at different laboratory angles as indi- 



cated in the figure. Full boxes with error bars are the experimental data taken from Ref. [25| and 



the results of QMD + SDM are denoted by solid histograms. The dashed histograms are results of 
another QMD + SDM obtained with only the resonance-like angular distribution Eq. (21), while 



the dot-dashed histograms with the direct-like angular distribution Eq. (|24|), respectively. 

FIG. 10. Neutron energy spectra for the reaction p(113 MeV)+ 56 Fe at different laboratory 
angles as indicated in the figure. The x-axis is plotted in a logarithmic scale in the left-hand-side, 
while the linear scale in the right-hand-side. The solid histograms are the results of QMD + SDM, 
the open circles with error bars denotes the experimental data taken from Ref. [^6[. The dashed 
histograms denote the results of NUCLEUS [§7) at the 150° laboratory angle. 
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FIG. 11. Same as Fig. M for the reaction p (256 MeV)+ 208 Pb. 



FIG. 12. Same as Fig. [H] for the reaction p (597 MeV)+ 56 Fe. 



FIG. 13. Same as Fig. 1C for the reaction p (800 MeV)+ 2U8 Pb. 



FIG. 14. Same as Fig. 10 for the reaction p (1500 MeV)+ 208 Pb. The full boxes with error bars 
are the experimental data taken from Ref. p3|. 



FIG. 15. Same as Fig. || for the reaction p (3000 MeV)+ 208 Pb. 
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TABLE I. Elastic cross section parameters. 






p-n 


others 


Ci (mb) 


28.0 


35.0 


C 2 (mb) 


27.0 


20.0 


C 3 (mb) 


12.34 


9.65 


C 4 (mb) 


10.0 


7.0 


TABLE II. Inelastic cross section parameters. 




fn cio 


<TN* 


a 


3.0 14.0 


23.0 


a 


0.9 -0.3 


1.5 


m (MeV) 


1188 1245 


1472 


r (MeV) 


99.02 120.0 


300.0 


TABLE III. Parameters in the width of resonances. 




M r (MeV) T r (MeV) 


r (MeV) 


A 


1232 110 


300 


N* 


1440 200 


523 


TABLE IV. Parameters in f R (s, cos 9). 


V^(GeV) 


y/s<2.U 2.14 <^/i< 2.4 


2.4 < 


ai(s) 


0.5 29.03 - 23.75s + 4.865s 2 


0.06 


03(5) 


0.0 -30.33 + 25.53s - 5.301s 2 


0.4 



36 



60 







i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 
(a) p-p Cross Sections 
total 




I I UlLJ L ■!■ ■!• ± it J I I I I I I I I I I I I I I I I I I 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

T L (GeV) 



60 



20 







i i i I i i i i I i i i i I i i i i I i i i i I i i i i 
(b) p-n Cross Sections 

total 



-O U WJ Ub UJuOU (H Off 

o o og o 




elastic 



inelastic 



I i i^kY I. .I. .!■ -i- •!' 'l I I I I I I I I I I I I I I I I I I 



r 



N 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

T L (GeV) 



20 



15 - 



-Q 



i i I i i i i I i i i i I i i i i I i i i i I i i i i 
NN -> NN* -> NArc -> NNrar ( o N - ) 

§ pn -> pp 7T° 7t" X 1 5 

4 pn -> pn 7r + 7r" x 1 2 / 5 

£ pp -> pp 7T + 7t" X 3 



10 



- o 



Q j^P I i 




i I i i i i I i i i i I i i i i I i i i i 



0.5 1.0 1.5 2.0 

T L (GeV) 



2.5 3.0 




cos(e) 




50 - 



O 

Q 

o 

DC 







CD 
O 

c 

CD 
CD 



-50 - 



-100 



— i 1 1 1 1 1 1 — 

Directed Transvererse Momentum 



„ -o- 



_ -E- 



_ — B 



_ _ E] 




RQMD 

— ■— Present QMD 

— g— Standard QMD N> 
- e - Standard QMD with Lorentz Contraction 

i I i I i I i 



o 



1 2 3 

E (GeV/ nucleon) 



(a) Density distribution 



40 Ca 



QMD 

HF 

Matter 



N i 



r (fm) 



15 



10 



Cl 







i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — r 

(b) Momentum distribution 

|~b— QMD 

HF 

Matter 



\ 



\ 



\ 



J I I I I I I I I I I L 




J 1 L 



0.0 0.5 1.0 1.5 2.0 

p(1/fm) 



10 4 i 



1 — I I I I 1 1 1 1 1 — I I I 1 1 1 1 1 1 — I I I I 1 1 1 1 

208 Pb(p,xn) E p = 1 500 MeV 



S Ishibashi et al. 
I - ^ QMD+SDM 
r — from QMD 
r ~~ - from SDM 




10 



J I I L 



i i n~i+ 



10° 



10 1 10 2 
E n (MeV) 



10' 



10 1 
10° 
10" 1 



O o 

o 
> 

w 



CO 

Q. 

D 10" 5 

CO 

LU 



10" 6 
10" 7 
10" 8 



n 1 i 1 r 
3.17 GeV 

- QMD+SDM - 
En'yo et al. 




E 45.0 

60.0° (x10 1 ) 
72.5° (x10 2 ) 
84.3° (x10 3 ) 

_J I I I L 



200 400 600 800 1000 1200 
E p (MeV) 



t | i | r 



1 i i i i r 



AI(p,X7i) 3.17 GeV 



QMD 

En'yo et al. 




45.0 

10 - 7 U 60.0° (x10 1 ) 
72.5° (x10 2 ) 
10 -8 L 84.3° (x10 3 ) 

_j I i I 



200 400 600 800 1000 1200 
E* (MeV) 



"1 1 — I — I I I I I I 1 1 — I — I I I I I 



56 



> 

CD 



j2 



10 £ 



10 1 , 



10° 
10" 1 



lu 10' 

a 

~B -, n -3 



10" 4 
10" 5 



Fe(p ; xn) E p = 113 



10 



-6 




QMD + SDM 
Meier et al. 
NUCEUS 



J I I I L 



10 c 



10 1 

E n (MeV) 




-| Q-6 I I I I I I I I I I I I I I I I 

20 40 60 80 100 120 
E n (MeV) 





1Q-9 I I I I I I I I I I I I I I 

100 200 300 

E n (MeV) 



t — i — i rrrra 




10° 10 1 10 2 10 3 

E n (MeV) 




Iq-8 I i LU i U i I l 

200 400 600 

E n (MeV) 



